Cooperative spontaneous emission in nonuniform media 
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Abstract. The subject of this paper is modification of cooperative spontaneous emission by a nonuniform 
medium, with nonuniform distributions of electromagnetic field. A brief analyzis is presented and it is 
postulated, that if spontaneous emission from an atom is strongly suppressed, cooperative emission with 
another atom may be a preferred emission channel and counteract the suppression. 
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1 Motivation 

It is known, that interaction of an atom with electromag- 
netic field is affected by tlie environment surrounding the 
atom, whether it is a resonant cavity pi or, generally, 
a nonuniform medium. In particular, modification of spon- 
taneous emission occurs, what has been first explained 
by Purcell [5]. This fact has found many practical ap- 
plications, e.g. following the suggestion of Yablonovitch, 
that a medium with a photonic band gap for all direc- 
tions of propagation, nowadays called a photonic crys- 
tal, could suppress spontaneous emission '3 . Another ef- 
fect concerning spontaneous emission, known as superra- 
diance, has been predicted by Dicke [3j and later confirmed 
in experiments: systems of many atoms can emit photons 
cooperatively, what results in radiation intensity propor- 
tional to N'^. Obviously, these effects combine, and in re- 
sult superradiance is modified by the environment as well. 
This paper is devoted to analysis of this modification in 
a nonuniform material, in which modification of sponta- 
neous emission is significantly different in various regions. 
Special attention is paid to the case of an atom is placed 
in a position in which spontaneous emission is strongly in- 
hibited. It is postulated, that cooperative emission might 
be a preferred decay channel in this case. 



2 Quantum mechanical description of 
spontaneous emission 

Theory of spontaneous emission is described in many books, 
e.g. [S], and papers, e.g. the one by Glauber and Lewen- 
stein [B], where a formulation proper for nonuniform di- 
electric medium is presented. The most important effects 
imposed on spontaneous emission by a material or a struc- 
ture can be described in a simple non-relativistic approach. 



The whole system is decomposed into two subsystems. 
The first one is the electromagnetic field, described by the 
hamiltonian 

iJem-^fiWA^A, (1) 
A 

where A runs through all the modes of the field, ujx is angu- 
lar frequency of the mode A and A^a is the photon number 
operator for this mode. Eigenstates of Nx form the Fock 
basis, which, for spontaneous emission, can be restricted 
to states with no photons |0) and states with one pho- 
ton I 1a). Modes of the structure can form a discrete or 
a continuous (or partly discrete, partly continuous) set, in 
this latter case summation should be replaced by integra- 
tion over proper quantities, represented by A - this is only 
a matter of notation, for simplicity notation for discrete A 
will be used here. 

The second part of the system is the atomic part, which 
is assumed to consist of a number of identical two-level 
atoms. For simplicity, any dependence of the environment 
on atomic states is neglected, assuming identical excita- 
tion energy Hf] and transition dipole moment d for each 
atom. State of A^at-atomic system can be defined as tensor 
product of individual atomic states 
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where |-l-„) and | — 
of the n-th atom, 
comes then 



^) denote the excited and ground state 
respectively. Atomic hamiltonian be- 
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To conveniently describe atomic transitions, one can 
define a relaxation operator for the n-th atom R„, corre- 
sponding to deexcitation of this atom; 
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Then, its hermitean conjugate Rjj corresponds to putting 
the n-th atom in the excited state. The interaction hamil- 
tonian allowing to describe spontaneous emission can be 
written in the form: 



^int = -fi^^gx^n{iRi ■ f\{r„)ax + h.c.) 



(5) 



A n=l 



where r„ is the position of the n-th atom, f\{r) is the elec- 
tric field distribution of mode A, abbreviation h.c. stands 
for hermitean conjugate of the previous term, and 

gx,n = 9\{rn) = Aici(r'„) ^.^ (6) 

The factor £dici('") describes relation between amplitudes 
of local and macroscopic electric field OIZ], which may 
be position dependent. a\ is the anihilation operator for 
mode A, which satisfies bosonic commutation rules with 



where C^^ stands for the inverse Laplace transform. The 
solution can be easily written as a function of time in 
Weisskopf-Wigner approximation, which is obtained by 
evaluating the second term in the nominator of C(s) in 
the limit s — > 0"*". This procedure leads to the well-known 
result 

|Co(t)|' = e-", (12) 
with the total decay rate 



t. 



its hermitean conjugate, the creation operator a^^. 
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+ E E C^,-A(t)c-'[(^-i)^+"^l* |1a) If) , (9) 

A *" 

where ^' runs through states with — 1 excited atoms. 
Initial conditions relevant for spontaneous emission are: 
Co(0) = 1 and the remaining amplitudes C^plx{Q) = 0. 
Projecting the Schrodinger equation onto individual states 
from the state vector one obtains a set of equations of mo- 
tion for amplitudes Cq and Gfi\, which are most elegantly 
treated by application of Laplace transformation, leading 
to a set of linear equations for their Laplace transforms. 
This way, one easily obtains a formal solution for Cq: 



(10) 



with 
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r = J2rx5{{2-ujx) 



(13) 



and rate of emission into mode A given by: 

2^ E 5A.n E 9x,n' {Ri ■ /A(r„)R„, • flirn,)) , (14) 



n—1 n' — l 



Let the system be initially in a state |0) Wo), where Wq) 
is a state with N excited atoms. Because coupling of the 
atomic system with electromagnetic field is weak, it is rea- 
sonable to consider only the states, which couple directly 
to |0) Wo), i-e. these, which combine into Hint |0) Wo)j which 
are states with one photon and — 1 excited atoms. All 
the other states can be produced by coupling of higher 
order, several magnitudes smaller, and will be neglected. 
This approximation allows to define the Hilbert space for 
the considered problem and conduct derivation of spon- 
taneous emission rate the same way, as for a relatively 
simple single atomic case [5]. State vector of the system 
can be defined as: 



with (•) denoting mean value in the state Wo)- This is 
a standard derivation, but its reconstruction helps to pre- 
cisely define the considered problem. 



3 Cooperative spontaneous emission 

The result ([H)) for the emission rate explains existence of 
the most important phenomena. Except for description of 
the simplest case of only one atom, it allows to observe 
the term proportional to N'^ and it reveals dependence 
of the decay on distribution of electric field. Thus, it is 
sufficient for an analyzis of spontaneous emission decay 
in structures, where distributions fx depend strongly on 
position, what is the case of interest in this paper. 
Let us split the expression (jl4|) into two terms: 
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(15) 



the first term describing contributions of individual atoms: 

r|^V) = 27r.gl„ (R,t • /A(r„)R„ • /^(r„)) (16) 

and the second one describing contributions of atomic 
pairs: 

27Tgx,n9x,n' ' /A(r-„)R„, • /l(r„0) • (17) 

These expressions allow to study how field distributions 
affect the decay rate. Contributions of individual atoms 
are modified accordingly to \f x{'r'7i)\^ ■ In particular, if no 
modes with angular frequency f2 are present in the struc- 
ture, emission at this frequency is completely stopped. 
This is the idea which has been presented by Yablonovitch 
[3]. However, this is impossible in real structures, which 
are of finite size and there always exists a set of modes 
with a given frequency. It is possible though, that a struc- 
ture very strongly extinguishes field of a mode in a certain 
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region. Therefore, strong inhibition of spontaneous emis- 
sion into a mode can occur, but because the structure 
adequately forms the field distribution, and not because 
this distribution is not present at all. 

This fact may have considerable consequences for col- 
lective emission by atomic pairs, which is described by the 
term {n,n'). If both atoms of a pair are placed nearby, 
they interact with modal fields having approximately the 
same intensities at their positions, thus, modification of 
emission is similar to single-atomic case and, apart from 
superradiant emission ^4, , no new effect is observed. How- 
ever, if the atoms are placed in positions, at which inten- 
sity of a mode |/a| differs significantly (this might occur 
e.g. in a finite one-dimensional photonic crystal for a mode 
from a photonic band gap, for which field distribution is 
strongly suppressed in the middle of the structure, but 
outside forms a partly standing wave with much larger 
amplitude W), rate of cooperative emission by the pair 
into the mode is proportional to product of these intensi- 
ties, significantly different than their second powers. Thus, 
for atoms n-th and n'-th, if |/a(t"„)| <C |/A(»"n')|j emission 

rate from the n-th atom into the mode r^^\n) is propor- 
tional to \f\{rn)\^ and can be significantly suppressed, 
but the factor related to electromagnetic field in the rate 

(2) 

of collective emission by the pair (n, n') is much big- 
ger |/A(i"n)| |/A(^n')l- This meaus, that electromagnetic 
field strongly privileges this emission channel and makes 
emission more probable. 

For two atoms to emit a photon cooperatively a corre- 
lation (at least partial entanglement) between their states 
is necessary. This correlation may appear beacause of ra- 
diative processes [9l[T^, thus cooperative emission may 
occur even if atoms are pumped incoherently. Another re- 
striction on cooperative emission is the existence of a max- 
imal distance, at which the atom can emit in cooperation: 
radiation emited by one of the atoms must reach the sec- 
ond atom before the act of emission is finished. This con- 
dition for atoms in vacuum has been discussed by Arecchi 
and Courtens [TT] . They have argued, that if the lifetime 
of the excited state is Tox, then the critical distance is sim- 
ply Lc = cTox and have provided a formula, which can be 
written in the following form: 

where Xq is the wavelength of the emitted radiation, J\f — 
concentration of radiating atoms and Fi - emission rate 
for a single atom. For optical transitions in gases Lc is 
of the order of 10 cm. Formula (flS)) can be used also for 
a rough estimation of Lc in a structure. Let us assume 
Ao « 1 fJ.rn, concentration J\f « 10^^ cm~^, and a value 
typical for optical transitions Fi « 10^ s~^. Finally, to 
conform to the idea of the condition, instead of c one 
should use Vg, group velocity of light in the structure. 
This leads to the result Lc ~ 3-y/ug/c x lO"'^ m. Even if 
the group velocity in the structure dropped to a few m/s, 
i.e. Vg/c « 10^*, it would lead to Lc ~ 30 ^m. More real- 
istic values can be found in literature, where observations 



of slow light with Vg/c w 10^^ are reported [T^HTl] and 
this ratio can be used as a representative value. Then, 
one obtains ic ~ 3 cm, a distance longer, than dimen- 
sions of optoelectronics devices. Therefore, the described 
phenomenon of collective emission into a mode, for which 
emission from a single atom would be strongly suppresed, 
could become a dominant process and might be observed 
in structures, where strong inhibition of spontaneous emis- 
sion occurs. This fact would have practical consequences, 
because it would counteract stopping the emission or en- 
hancing the lifetime of atomic excited state by a structure 
strongly modifying distributions of electromagnetic field. 

Physical systems are more complicated than just iso- 
lated ensembles of atoms and, in general, interactions of 
different kind than direct radiative coupling between the 
atoms are present. These interactions can also cause col- 
lective states to appear and contribute to superradiant 
emission this way. For example, recently existence of su- 
perradiance in a system with plasmonic coupling has been 
discussed 15i . A similar effect of structure nonuniformity 
should be observable in these systems as well. 
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